Low temperature hep to monoclinic structural transition in solid C70 : Ephemeral 

nature of the intermediate phase 
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We follow the structural transformation in solid C70 from the high temperature hep to a low 
temperature monoclinic phase using detailed x-ray diffraction studies at controlled cooling-rates from 
0.0033 to 0.42 K/min. Rapid cooling of the sample gives the signature of a two-step transformation 
which disappears when the system is transformed quasi-statically. These experimental results can 
be rationalized using a simple mean field, Langevin dynamical theory using a free energy functional 
with minima corresponding to the parent and two competing product phases such that one of these 
product phases remains metastable throughout. The implication of our results on the existence of 
the intermediate phase in the sequence, hep-dhep- monoclinic, of structural transitions in solid C70 
with the lowering of temperature is discussed. 
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Following the pioneering work of Heiney et ai, ]jj on 
relating the structural transition at low temperature (260 
K) in solid Ceo to the development of orientational corre- 
lations between rotating Ceo molecules, there have been 
extensive studies on fullerene molecular crystals, both 
solid Ceo and C70. However, unlike in solid Cqq, the 
low temperature structural transitions in the next higher 
fullerene C70 are not unambiguously established. The 
situation is exacerbated by the complicated orientational 
dynamics of the lower symmetry (Dsn) molecule, asso- 
ciated sluggishness of the transition and complications 
due to intercalated solvent effects etc. Solid C70 is ob- 
served to crystallize in both fee and ideal hep structures, 
corresponding to close packing of rapidly rotating spher- 
ically averaged molecules. The fee solid C70 is first seen 
to transform at «280 K to a rhombohedral structure 
and subsequently to a monoclinic structure at 200 K || 
with the lowering of temperature; both these transitions 
are associated with considerable hysteresis of « 50 K. 
A widely accepted viewpoint on the nature of structural 
transitions in hep solid C70 involves two transitions : One 
from the ideal hep (c/a w 1.63) to the deformed hep (c/a 
w 1.82) at 337 K (dhep) and a second transition to the 
monoclinic structure at 276 K. ||-|5| Data from differ- 
ential thermal analysis (DTA) differential scanning 
calorimetry (DSC) and dilatometry measurement ||] 
also lend support to this scenario. The occurence of a 
two-step transition is rationalized in terms of successive 
freezing of molecular orientations, first along the short 
C2 axis and then along the long C5 axis; the structural 
transitions being produced by the coupling of these ori- 
entational degrees of freedom to the lattice strain. 

Several experiments, however, see a single-step transi- 
tion from hep to monoclinic phase, and there are others 
Q that have reported three structural transitions with 
the variation of temperature. For example, Dennis et al. 
pot observed a single-step orientational ordering of C70 
molecules, over the temperature range of 270 to 160 K 
in their muon spin resonance (/iSR) experiments. In our 
earlier work jllj , using high quality samples prepared by 
the solution route we observed (1) the equilibrium room 



temperature structure of solid C70 to be ideal hep and 
(2) a broad (270 to 170 K) single-step transition from 
the ideal hep to a monoclinic structure implying a si- 
multaneous freezing of all the orientational degrees of 
freedom. It is of interest to know if the differences in 
the structural sequence are related to the different (of- 
ten uncontrolled) cooling-rates employed in the earlier 
experiments. If an equilibrium intermediate phase actu- 
ally exists then it would be seen during controlled slow 
cooling of the system from the ideal hep structure. 

With a view to investigate this, x-ray diffraction 
(XRD) experiments at various controlled cooling-rates 
(varying over two orders of magnitude) have been car- 
ried out, and the evolution of the monoclinic phase at 
low temperature followed. Contrary to our initial ex- 
pectations, though the transformation proceeds through 
a single-step (starting at «270 K) when a slow cooling- 
rate (0.0033 K/min.) is used, it has a two-step sequence 
on cooling rapidly (0.42 K/min.). Our experimental re- 
sults (Fig. 2.) strongly suggests that the intermediate 
phase reported earlier is likely to be an artifact of 
the kinetics of the transformation process and not a true 
equilibrium phase. We show that qualitative aspects of 
our data for the temparature evolution of the fraction of 
the retained hep phase during the hep to monoclinic tran- 
sition can be understood within a simple phcnomcnolog- 
ical kinetic theory employing a general free energy func- 
tional with minima corresponding to the parent and two 
product phases such that one of the product phases is 
always metastable. 

Chromatographically separated C70 powder was 
heated in vacuum at 250° C for 24 hours to get rid of 
trapped solvents. The x-ray diffraction (XRD) study at 
low temperatures was carried out using a Siemens D-500 
powder diffractometer with a continuous flow cryostat. 
The sample holder was a (911) cut silicon wafer, JL2) at- 
tached to the cold tip of the cryostat using a thin layer of 
GE vernish. The sample displacement errors due to ther- 
mal contraction were corrected by simultaneously using 
a small quantity of NIST 640a silicon powder as internal 
d-spacing standard. The XRD patterns were recorded at 
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10 K steps over the temperature range of 300 to 77 K. 
The sample was cooled at controlled rates. Experiments 
were carried out at four different cooling-rates varying 
over two orders of magnitude, viz. 0.0033 K/min., 0.22 
K/min., 0.31 K/min., and 0.42 K/min. Before commenc- 
ing a fresh cooling cycle we carried out a detailed XRD 
run at room temperature to ensure that the sample had 
regained its original hep structure. 
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Fig. 1 XRD patterns showing the temperature evolu- 
tion of the monoclinic phase from the parent hep at different 
cooling rates, (a) 0.0033 K/min.; the lowest cooling-rate em- 
ployed in this experiment (b) 0.22 K/min. (c) 0.31 K/min. 
and (d) 0.42 K/min. 



change commences at ?s270 K, and by sal70 K the trans- 
formation to the monoclinic phase is complete. While 
apriori, such a broad transition could be attributed to 
impurity effects (presence of intercalants) , the observa- 
tion of a «80 K co-existence region even in the sublimed 
sample Q suggests that this feature is intrinsic. It is seen 
from Fig. 2 that the width of the transition increases with 
the increase in cooling rate (e.g., 0.22 K/min.) implying 
that the width is related to the sluggishness of transition 
in solid C70. The kinetic model proposed in this pa- 
per (see below) accounts for the broad hep to monoclinic 
transition solid C70 (c/. Fig. 2 and Fig. 3). 

The interesting aspect of Fig. 2 is that with increas- 
ing the cooling-rate further, one obtains a step in the 
transformation region, which can be seen as a distinct 
plateau for the cooling-rate of 0.42 K/min, the fastest 
controlled cooling rate possible with our present experi- 
mental setup. 
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As shown in Fig.l (a), the room temperature XRD 
pattern can be indexed to the ideal hep with lattice pa- 
rameters a = 10.593 A , c = 17.262 A , and space group 
P6j,lmmc 13 1. The low temperature pattern can be in- 



dexed to a monoclinic structure, space group P2\/m 1 
with lattice parameters a= 10.99 A , b— 16.16 A , c=9.85 
A ,13= 107.75°. At a slow cooling rate of 0.0033 K/min, 
it is seen that the hep phase transforms to the monoclinic 
structure without the intervening dhep phase, consistent 
with the results reported earlier pj. At higher cooling 
rates of 0.22, 0.31 and 0.42 K/min, XRD patterns in the 
29 range of 16° and 18°, corresponding to the most in- 
tense peaks of the hep (110) and the monoclinic (211) 
phases |0], have been measured, and these results are 
shown in Figs.l (b)-(d). 

The change in the fraction of the hep phase with tem- 
perature has been calculated from these patterns as the 
area ratio of the 100% hep (110) reflection to the sum 
of the areas of the (110) reflection of the hep and (211) 
reflection of the monoclinic phase, and these results are 
plotted in Fig. 2. It is seen that at a cooling-rate 0.0033 
K/min. the hep fraction decreases in a single-step, within 
a temperature range of about 100 K. The structural 



Fig. 2 Change in the fraction of the retained parent 
hep phase with temperature at different cooling-rates as in- 
dicated. 

We argue below that this plateau is a signature of the 
formation of an intermediate dhep phase during the hep 
to monoclinic transition. The dhep phase has been re- 
ported § to have lattice parameters of a = 10.11 A , c = 
18.50 A . The (103) reflection of this structure occurs at 
17.59, very close to the (211) line of the monocline phase 
occuring at 17.61, and has an intensity that is only w 30 
% intensity of the strongest (110) reflection ||. Thus, the 
formation of the dhep phase will result in a decrease (in- 
crease) in the estimated monoclinic (hep) fraction, lead- 
ing to the plateau. The interesting aspect of Fig. 2. is 
that this plateau, corresponding to the dhep phase, is 
seen during fast cool. If the intermediate dhep phase 
were an equilibrium phase, we should have observed it in 
the slow cool experiments. This suggests that interme- 
diate phase is metastable, and the feature seen in Fig. 2 
is a reflection of kinetic effects - a model for which is 
developed as follows. 
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Being a large and somewhat asymetric molecule C70, 
possesses a rich energy landscape with many closely com- 
peting local minima with the possiblity of the system get- 
ting trapped in metastable phases influencing the equilib- 
rium phase transitions in this system. A Landau theory 
for the hep - monoclinic transition in C70 can, in prin- 
ciple, be constructed along the lines worked out in Refs 
for the fee - rhombohedral transition. In this case, 
the free energy would have three nonequivalent minima 
corresponding to the three (hep, dhep and monoclinic) 
phases described above and is expected to have either 
one (hcp-monoclinic) or a sequence of two (hep-dhep- 
monoclinic) first order transitions depending on the rel- 
ative stability of the dhep and the monoclinic phases. 
The question of relative stability is, unfortunately, tricky 
since it depends on the actual magnitudes of the var- 
ious coefficients in the Landau expansion. The uncer- 
tainty in the values of these coefficients, which have to 
be obtained experimentally or from accurate microscopic 
calculations, makes it difficult in practice to settle this 
question convincingly. 
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Fig. 3 The time and noise averaged order parameter I-77 
plotted as a function of the dimensionless temperature r for 
various cooling rates obtained by holding the system for th 
steps as t is reduced from r = 1 in steps of .01. The order 
parameter r\ is obtained by numerically solving the stochastic 
differential differential equation Eq.3. (inset) Plot of the free 
energy f(rj) as a function of the order parameter 77 for r = .1, 
and — .1. Note that the intermediate minimum at 77 = 0.4 is 
always metastable. 

We therefore choose a simple, phenomenological ap- 
proach to understand the qualitative features of the ki- 
netics of the hcp-monoclinic transition as depicted in Fig. 
2. without trying to derive them with any amount of nu- 
merical accuracy. Consider a general free energy f(rj) 
with three minima obtained by adding inverted Gaus- 
sians at r\=r\\ and 772 with widths ad and a2 representing 
the product phases, to a parabolic potential centred at 
rj = 0, representing the parent phase. 



/ = r/ 2 {l - (1 - r)((5e - ° l(lJ -'» l)a + eT^"-"*) 2 )} (1) 

The depth of these minima relative to the one at 77 = is 
determined by the dimensionless "relative temperature" 
t (x (T — T c )/T c (defined for r < 1 with T c the equilib- 
rium transition temperature) and 8 controls the relative 
stability of the products for 771 and 772; we have chosen 8 
to be such that the minimum at 771 is always metastable 
throughout the range of r considered. At t = we have 
a single first order transition from the parent to the prod- 
uct corresponding to 77 — r] 2 - Using the (arbitrary) set 
of parameters 771 = 0.4, a% = 200, 772 = 0.8, a 2 = 10 and 
8 = 0.78 we have plotted the free energy f(rf) for r = .1, 
and —.1 in Fig. 3 (inset). Note that the parameter 8 has 
been so chosen that the intermediate minimum is always 
slightly metastable compared to the other two minima. 
While our order parameter 77 has no direct connection to 
actual order parameters for orientational transitions of 
C70 it may be thought of as a "distance" along a "mini- 
mum free energy" path connecting the parent minima to 
the two nonequivalent product minima in the true high 
dimensional order parameter space. 

We now investigate the dynamics of this transition 
and especially the effect of the metastable minimum at 
?/ = 7 n. The dynamics is governed by a Langevin equa- 
tion [|15| for the order parameter 77 (we neglect all spatial 
fluctuations) evolving in the free energy landscape f(rj) 
and is given by, 

dr 1 /dt = -Tdf{r,)/dr 1 + m- (2) 

The friction coefficient T is related to the variance of 
the Gaussian white noise 7 by the usual fluctuation- 
dissipation relation 7 =< >= 2k B TTS(t - t'), 
where ks is the Boltzmann constant and T the actual 
temperature. The value of T can, in principle, be ob- 
tained from the libron lincshape obtained from Raman 
fl6|| and neutron jl7) scattering experiments. The angu- 
lar brackets signifies an average over various realizations 
of the random noise. 
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Fig. 4 The normalized probability distribution P(rj, r(t)) 
of the order parameter r\ as a function of 77 as r is reduced 
from t = 1 with a holding time of th = 100 steps (curve 
number 3 in Fig. 3) 

We simulate Eq.|| using a straightforward Euler dis- 
cretization scheme with a timestep of At = .01. The 
dimensionless temperature r is reduced from 1 in steps 
of .01. At each temperature step the system was held for 
th steps (a larger value of th implies slower cooling) . The 
order parameter rj averaged over th steps and over vari- 
ous realizations of the random noise (typically 300 — 700) 
is plotted as a function of r in Fig. 3 (for the above 
set of parameters for f(rj) and 7 = 1.8). The resem- 
blance with the cooling curves of Fig. 2 is remarkable. 
Both extremely fast and slow cooling obtains a single step 
transformation — the transformation becoming sharper 
as the cooling rate is reduced. For intermediate values of 
the cooling rate a plateau develops. The plateau results 
when the system is trapped, temporarily, in the interme- 
diate phase as it evolves from the parent to the product 
structure. This can be illustrated as follows. 

In Fig. 4 we have plotted the probability distributions 
P(rj,T(t)) as a function of r(t) for a collection of 500 in- 
dependent systems cooled from r — 1 to t — —4.0 with 
th = lOOAi. Starting from a delta function at the origin 
at t = the P(j], t) spreads out with time and eventu- 
ally evolves into a strong peak centered at the equilib- 
rium value of r\. However, we observe that the system 
spends a considerable amount of time in the secondary 
metastable minimum near rji which roughly corresponds 
to the plateau observed in Fig. 3 (curve 3). In other 
words, at slow cooling rates the system gets enough time 
to relax to the global minimum. For very fast quenches, 
on the other hand, by the time the system relaxes from 
the parent to the metastable intermediate phase, r de- 
creases to the extent that the barriers are no longer suf- 
ficient to trap it there. Similar qualitative behaviour is 
seen over a broad range of values for our parameters as 
long as the intermediate minimum lies sufficiently close to 
the true global minimum and it is sharp enough to trap 



the system effectively. This behaviour is quite general 
and is analogous to that in ferroelectrics where a similar 
kinetic stabilization of the metastable paraelectric state 
produces double hysteresis loops [[L8|| . 

To paraphrase, the low temperature structural trans- 
formation of solid C70 exhibits interesting kinetic ef- 
fects, as evidenced by XRD experiments carried out at 
various cooling rates. At extremely slow cooling rates 
(0.0033K/min), rarely employed in the studies so far, the 
transformation to the monoclinic phase proceeds through 
a single step. This implies that the equilibrium trans- 
formation corresponds to the simultaneous freezing of all 
orienational degrees of freedom, and recently |r]] a struc- 
tural model for the low temperature monoclinic structure 
has been proposed. At high cooling rate the transforma- 
tion from ideal hep to monoclinic phase proceeds through 
an intermediate phase. With support from a theoretical 
model, we show that such a situation can arise if the in- 
termediate phase is metastable. In the light of present 
investigations theoretical calculations on the stabilty of 
dhep phase need to be re-examined. 
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